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Abstract 

"q^" ■ We study angularly excited as well as interacting non-topological solitons, so-called Cu- 

bans and their gravitating counterparts, so-called boson stars in 3 + 1 dimensions. Q-balls 
. and boson stars carry a non- vanishing Noether charge and arise as solutions of complex scalar 

field models in a flat space-time background and coupled minimally to gravity respectively. 
We present examples of interacting Q-balls that arise due to angular excitations, which are 



closely related to the spherical harmonics. We also construct explicit examples of rotating 
boson stars that interact with non-rotating boson stars. We observe that rotating boson 
stars tend to absorb the non-rotating ones for increasing, but reasonably small gravitational 
coupling. This is a new phenomenon as compared to the flat space-time limit and is related 
to the negative contribution of the rotation term to the energy density of the solutions. 
1} | In addition, our results indicate that a system of a rotating and non-rotating boson star 



can become unstable if the direct interaction term in the potential is large enough. This 
instability is related to the appearance of ergoregions. 
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^9 ' 1 Introduction 

On 

co 

. Solitons play an important role in many areas of physics. As classical solutions of non-linear 

field theories, they are localised structures with finite energy, which are globally regular. In 
general, one can distinguish topological and non-topological solitons. While topological solitons 
. ^ | [lj possess a conserved quantity, the topological charge, that stems (in most cases) from the 

spontaneous symmetry breaking of the theory, non-topological solitons [2l [3] have a conserved 
S ■ Noether charge that results frora a s y mmet ry of the Lagrangia, The standard exarapfe of 

non-topological solitons are Q-balls 0j, which are solutions of theories with self-interacting 
complex scalar fields. These objects are stationary with an explicitely time-dependent phase. 
The conserved Noether charge Q is then related to the global phase invariance of the theory and 
is a function of the frequency. Q can e.g. be interpreted as particle number [2]. 

While in standard scalar field theories, it was shown that a non-renormalisable $ 6 -potential 
is necessary [3], super symmetric extensions of the Standard Model (SM) also possess Q-ball 
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solutions [6]. In the latter case, several scalar fields interact via complicated potentials. It was 
shown that cubic interaction terms that result from Yukawa couplings in the superpotential and 
supersymmetry breaking terms lead to the existence of Q-balls with non-vanishing baryon or 
lepton number or electric charge. These supersymmetric Q-balls have been considered recently 
as possible candidates for baryonic dark matter j7] and their astrophysical implications have 
been discussed [8j. In [9], these objects have been constructed numerically using the exact form 
of the supersymmetric potential. 

Q-ball solutions in 3 + 1 dimensions have been studied in detail in [5] I10|, fTTj . It was realised 
that next to non-spinning Q-balls, which are spherically symmetric, spinning solutions exist. 
These are axially symmetric with energy density of toroidal shape and angular momentum 
J = kQ, where Q is the Noether charge of the solution and k £ Z corresponds to the winding 
around the z-axis. Approximated solutions of the non-linear partial differential equations were 
constructed in [5] by means of a truncated series in the spherical harmonics to describe the 
angular part of the solutions. The full partial differential equation was solved numerically in 
|10l llll [T2] (for a short review see [E]). It was also realised in [5] that in each fc-sector, parity- 
even (P = +1) and parity-odd (P = —1) solutions exist. Parity-even and parity-odd refers 
to the fact that the solution is symmetric and anti-symmetric, respectively with respect to a 
reflection through the x-y-plane, i.e. under 9 — > n — 9. 

These two types of solutions are closely related to the fact that the angular part of the 
solutions constructed in [5j [TU1 [TT] is connected to the spherical harmonic Yq(9,ip) for the 
spherically symmetric Q-ball, to the spherical harmonic Y^(9,ip) for the spinning parity even 
(P = +1) solution and to the spherical harmonic Y^(9, ip) for the parity odd (P = —1) solution, 
respectively. Radially excited solutions of the spherically symmetric, non-spinning solution were 
also obtained. These solutions are still spherically symmetric but the scalar field develops one 
or several nodes for r s]0, oo[. In relation to the apparent connection of the angular part of the 
known solutions to the spherical harmonics, "^-angular excitations" of the Q-balls corresponding 
to the spherical harmonics Y^{9, ip), —l<k<l have been constructed explicitely for some values 
of k and I in [12]. These excited solutions could play a role in the formation of Q-balls in the 
early universe since it is believed that Q-balls forming due to condensate fragmentation at the 
end of inflation first appear in an excited state and only then settle down to the ground state [H] . 
The fact that these newly formed Q-balls are excited, i.e. in general not spherically symmetric 
could, on the other hand, be a source of gravitational waves |15j . 

The interaction of two Q-balls has also been studied in [12J. It was found that the lower 
bound on the frequencies u>i, i = 1,2 is increasing for increasing interaction coupling. Explicit 
examples of a rotating Q-ball interacting with a non-rotating Q-ball have been presented. 
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Complex scalar field models coupled to gravity exhibit a new type of solution, so-called 
"boson stars" [T6J [TTl [18] . In pUHH] boson stars have been considered that have flat space-time 
limits in the form of Q-balls. 

In this paper, we study interacting boson stars with flat limit solutions in the form of the 
interacting Q-balls that have been studied in [12]. The boson stars are interacting via a potential 
term, but of course also through gravity. 

The paper is organised as follows: in Section 2, we give the model, Ansatz and boundary con- 
ditions. In Section 3 and 4, we discuss our results for Q-balls and their gravitating counterparts, 
boson stars, respectively. Section 5 contains our conclusions. 

2 The model 

In the following, we study a scalar field model coupled minimally to gravity in 3 + 1 dimensions 
describing two interacting boson stars. The action S reads: 

where R is the Ricci scalar, G denotes Newton's constant and C m denotes the matter Lagrangian: 

C m = -\d^ x d^\ - ^d^ 2 d^* 2 - $ 2 ) (2) 

where both $i and $2 are complex scalar fields and we choose as signature of the metric 
( — h ++). The potential reads: 

2 

v(*!, $ 2 ) = Yl N $ *l 6 - + 7il*il 2 ) + A l $ il 2 l$2| 2 (3) 

1=1 

where Ki, Pi, j{, i = 1,2 are the standard potential parameters for each boson star, while A 
denotes the interaction parameter. The masses of the two bosonic scalar fields are then given 
by (m^) 2 = 7i , i = 1,2. 

Along with [101 EH E2] , we choose in the following 

Kt = 1 , A = 2 , 7l = 1.1 , i = 1,2 . (4) 

In [5] it was argued that a <I> 6 -potential is necessary in order to have classical Q-ball solutions. 
This is still necessary for the model we have defined here, since we want $1 = and $2 = to 
be a local minimum of the potential. A pure $ 4 -potential which is bounded from below wouldn't 
fulfill these criteria. 
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The matter Lagrangian C m ([2]) is invariant under the two independent global U(l) transfor- 
mations 

$1 -» $ie ixi , $ 2 -» ^2e* X2 • (5) 

As such the total conserved Noether current j^ tot y \i = 0, 1,2, 3, associated to these symmetries 
is just the sum of the two individually conserved currents j± and j/j 1 with 

j(tot) = ii + # = - - i md^2 ~ *2d*** 2 ) ■ (6) 

with ft Jlt = 0, j$ ;M = and jg ot) ;/i = 0. 

The total Noether charge Q(tot) °f the system is then the sum of the two individual Noether 
charges Q\ and Q 2 : 

Q(tot) = Qi + Q2 = - J fl<$X - j j° 2 d 3 X (7) 
Finally, the energy-momentum tensor reads: 

2 

V = iWiWi + du^d^*) - g^C (8) 
i=i 

2.1 Ansatz and Equations 

For the metric the Ansatz in Lewis-Papapetrou form reads |10j : 



ds 2 = -fdt 2 + - (g{dr 2 + r 2 d6 2 ) + r 2 sin 2 6{d<p + ™dt) 2 ) (9) 

where the metric function f, I, g and m are functions of r and 9 only. For the scalar fields, the 
Ansatz reads: 

$i(t, r, 0, v?) = e^+^^r, 0) , i = 1, 2 (10) 

where the and the ki are constants. Since we require <&i(<p) = &i(<p + 27r) , i = 1,2, we have 
fcj 6 Z. The mass M and total angular momentum J of the solution can be read off from the 
asymptotic behaviour of the metric functions [lOj: 



M = — — lim r 2 d r f , J = —— lim r 2 m . (11) 

2G r^oo 2G r^oo 

The total angular momentum J = J\ + J2 and the Noether charges Qi and Q 2 of the two boson 
stars are related by J = k\Q\ + k 2 Q 2 . Boson stars with fcj = have thus vanishing angular 
momentum. Equally, interacting boson stars with k\ = — k 2 and Qi = Q 2 have vanishing angular 
momentum. 

The coupled system of partial differential equations is then given by the Einstein equations 

Gfiv = SttGT^ (12) 



with T^y given by ([8]) and the Klein-Gordon equations 

( D + W)^ = ° ' Z = 1 ' 2, (13) 
Details about these equations for one complex scalar field can e.g. be found in [TU] , 

2.2 Boundary conditions 

We require the solutions to be regular at the origin. The appropriate boundary conditions read: 

<9 r /| r=0 = 0, d r l\ r=0 = 0, g\ r=0 = 1 , m\ r=0 = , 4>i\ r =o = , i = l,2. (14) 

for solutions with ki ^ 0, while for ki = solutions, we have 

d r f\ r=0 = , <9 r /| r=0 = 0, g| r =o = 1 , m| r=0 = , <9 r 0i| r =o = , £ = 1,2. (15) 

The boundary conditions at infinity result from the requirement of asymptotic flatness and finite 
energy solutions: 

f\r— »oo — 1 j l\r— >oo — 1 i slr^oo — 1 j — >oo — i 0i|r^oo — , i — 1,2. (16) 

For = the regularity of the solutions on the z-axis requires: 

def\e=o = , d e l\g=o = , g| e=0 = 1 , <%m|e =0 = , <t>i\e=o = , £ = 1,2, (17) 

for k{ ^ solutions, while for fcj = solutions, we have 

«9e/|e=o = 0, d e l\ e=0 = , c/| e=0 = 1 , <%m| e=0 = , 9 e 0i| e=o = , £ = 1,2, (18) 

The conditions at 9 = ir/2 are either given by 

<%/|e =7r /2 = , del\ g =n/2 = , d e g\ e=n / 2 = , d e m\ e=n/2 = , <%^| e=7r/2 = , £ = 1,2 

(19) 

for even parity solutions, while for odd parity solutions the conditions for the scalar field functions 
read: 4>i\ e=7T / 2 = 0, £ = 1,2. 

2.3 Numerical procedure 

In the following, we will solve the system of partial differential equations (I12p and (|13p subject 
to the appropriate boundary conditions given in Section 2.2. This has been done using the PDE 
solver FIDISOL |19j . We have mapped the infinite interval of the r coordinate [0 : oo] to the 
finite compact interval [0 : 1] using the new coordinate z := r/(r + 1). We have typically used 
grid sizes of 150 points in r-direction and 70 points in 9 direction. The solutions pesented here 
have relative errors of 10 -3 or smaller. 
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Table 1: Boundary conditions and parity P for different choices of I and k 



3 Q-balls 

In this section, we study solutions in a flat space-time background, i.e. we choose G = 0. The 
Q-balls then interact only via the potential for A 7^ 0. 

3.1 Single Q-balls 

In order to be able to understand the structure of a system of two interacting Q-balls and their 
gravitating counterparts, we briefly reconsider the one Q-ball system. This has been studied in 
great detail in [5j [101 EJ 12] and arises from our model for A = and 02 = 0. In this section, 
we would like to emphasize some of the properties which are important to understand the 
gravitating case. In the present section on single Q-balls, we will omit the index 1, respectively 
2 for all quantities. As was noticed in [12] - using spherical coordinates (r, 9, 93) and a standard 
separation of variables - the solutions to the linearized scalar field equation are given by: 

( f>(r,e, i p)oc Jl+1/ ^ ujr) Y l k (e,<p) (20) 

where J denotes the Bessel function and Y, are the standard spherical harmonics with I integer 
and —l<k<l. As was shown in [12] - at least for the first few values of k and I - solutions 
to the full field equations exist in correspondence to the symmetries of the spherical harmonics. 
Explicit examples of the angular excited solutions have been presented for specific values of k 
and /. In the following, these solutions will be labelled by the two quantum numbers k and I, 
i.e. if we use we refer to the solution of the full equation which corresponds to the spherical 
harmonic . 

As already stated in Section 2.2, the boundary conditions depend on /, k. Here, we want 
to state the boundary conditions explicitely together with the parity P of the solution under 
— > tv — 9 for the first few values of k and I (see Table 1). Let us stress that all these conditions 
are compatible with the trivial solution = 0. Therefore, the numerical construction of non- 
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vanishing solutions turns out to be a non trivial task. Solutions for the choices of I and k given 
in Table 1 have been presented in [12J . 

3.2 Interacting Q-balls 

Interacting Q-balls were studied in detail in [12] . The interaction is characterized by the coupling 
constant A. The two Q-balls decouple in the limit A = 0. Here, we want to present the energy 
density of some of the solutions in order to illustrate the influence of the interaction term. 




Figure 1: The energy density of two interacting Q-balls with l\ = 0, k\ = (spherically sym- 
metric, non-rotating) and I2 = 1, ki = 1 (axially symmetric, rotating, parity even). Here A = 1 
and u\ = UJ2 = 0.8. Note that we use cylindrical coordinates z = rcosO and p = rsinO. 

In FigfT] we present the energy density for a spherically symmetric, non-rotating Q-ball 
interacting with an axially symmetric, rotating and parity even Q-ball for A = 1 and u\ = U2 = 
0.8. Note that we use cylindrical coordinates here with z = rcosO and p = rsin.6. The local 
maximum appears at z = 0, p ~ 2.2 and is connected to the rotating Q-ball. In Fig. [2] we give 
the energy density corresponding to a spherically symmetric, non-rotating Q-ball with l\ = 0, 
k% = interacting with an axially symmetric, non-rotating, i.e. angularly excited (J"t>all with 
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Figure 2: The energy density of two interacting Q-balls with l\ = 0, k% = (spherically sym- 
metric, non-rotating) and I2 = 1, &2 = (axially symmetric, non-rotating). Here A = 1 and 
uj\ = UJ2 = 0.8. Note that we use cylindrical coordinates z = rcos9 and /? = r sin 6*. 

Z2 = 1, = i s given. Again, we have chosen A = 1 and cj± = 002 = 0.8 for this case. The 
energy density has local maxima at three different values of z. The energy density thus consists 
of two tori at z > and z < and a deformed spherical ball at z = 0. 

4 Boson stars 

Once gravity is included (G / 0), the different possible Q-ball solutions discussed in the previous 
section become deformed by gravity and are called boson stars. Since only k 2 appears in the 
equation for the Q-ball the sign of k does not matter and we can restrict ourselves to k > 
without loosing generality in the flat space-time background case. This changes if we consider 
solutions in curved space-time, i.e. for G > 0. The reason is that the energy-momentum tensor 
involves terms that are linear in k. Reversing the sign of k the solutions will have the same 

1 Note that we plot the solution only for z > 0, but that we can continue it to z < due to its symmetry. 
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energy, but opposite angular momentum. 

In the case of interacting boson stars, the situation is more involved: when two boson stars 
with k\ and k 2 of the same sign interact the angular momentum is non-vanishing. However, for 
k\ = —k 2 the function m(r,8) = and the total angular momentum is vanishing. 

In the presence of gravity, the boson stars interact by means of the direct coupling term in 
the potential for A / but also through gravity. To understand the pattern of solutions we 
study the energy density in more detail. This reads 

r ° = vfa,<h) + e f g (dr^r + + -fn + (jpijj - j?) m (2i) 

The important point about this expression is that the term containing the "rotation function" 
m(r, 9) leads to a negative contribution to the energy density. Hence, rotation tends to decrease 
the energy of the solution. 

4.1 Numerical results 

In the following, we will abbreviate a = 8ttG and focus on the case u\ = uj%. We have fixed the 
parameters of the potential according to (J2j). 

4.1.1 h = 0, k x = 0, l 2 = 0, k 2 = 

This case describes two interacting non-rotating boson stars. The solutions are spherically 
symmetric. In particular, we find g(r) = 1 and m(r) = 0. The bound state of two spherically 
symmetric boson stars with a given value of a behaves roughly as a single boson star with 
a/2. The mass M and the values 0i(O), i = 1,2, /(0), 1(0) and Tq(0) slowly decrease when the 
parameter a increases. We found no evidence of a limiting behaviour, i.e. no critical a beyond 
which the solutions cease to exist. 

4.1.2 h = 0, h = 0, l 2 = 1, k 2 = 1 

This corresponds to the interaction of a non-rotating boson star and a rotating, axially symmetric 
and parity even solution. 

The contour plot of a typical solution of this type is presented in Fig. [3]fora = 0.2, A = 
and uj\ = a>2 = 0.8. Clearly, the metric functions show non-trivial behaviour now. 

Varying the gravitational constant a (with all other parameters fixed), our numerical results 
show that the scalar field function 4>i of the spherically symmetric boson star tends uniformly 
to zero and that only the spinning boson star survives for sufficiently large a. This result is 
consistent with the fact that rotation, i.e. non-vanishing m(r, 6) tends to decrease the energy 
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(see (|2ip ). The rotating boson star is therefore energetically favoured in comparison to the 
non-rotating case. This phenomenon is however strongly related to the gravitational interaction 
of the two boson stars. Indeed, the non-rotating spherically symmetric boson alone (for 02 = 0) 
exist for (arbitrarily) large values of a. This is illustrated in Fig. [3] for A = and uj\ = U2 = 0.8, 
where we plot the mass M, charge Q and </>i(0), which corresponds to the maximal value of 
the scalar field function <j>i{r) for non-rotating, spherically symmetric solutions The quantities 
corresponding to interacting boson stars are given in solid lines. <^i(0) and with this the field 
4>i(r) corresponding to the non-rotating boson star becomes identically zero for a = a cr ~ 0.25 
while the single boson star (dashed lines) exists for much larger values of a. To illustrate this 
phenomenon further, we also plot the energy density of these solutions for two different values 
of a in Fig. [5j The solution for a = 0.1 shows a quite large values for Tq at the origin p = 0, 
z = 0. Since a non- vanishing contribution to the energy density at the origin can only come 
from the spherically symmetric, non-rotating boson star, the contribution of this boson star is 
still quite strong here. When increasing a, the axially symmetric, rotating boson star tends to 
absorb the non-rotating one. This is clearly seen when comparing the a = 0.1 plot with that for 
a = 0.2. Here, the value of Tq at the origin has decreased strongly and the energy density tends 
to the shape of a torus signalling that the spherically symmetric contribution nearly vanishes. 

We observe the qualitatively same effect for A / 0. The critical value of the gravitational 
coupling a cr where <f>i becomes identically zero depends slightly on the coupling constant A. We 
find e.g. a cr (X = 0.5) rj 0.30, a cr (X = 0) ~ 0.25 and a cr (A = —0.5) rj 0.22, respectively. 

Since the appearance of ergoregions for globally regular solutions signals the existence of 
instabilities [20] we have studies these here as well. The ergoregions of single boson stars have 
been studied extensively in [TT]. Ergoregions exist if the goo component of the metric becomes 
positive, i.e. for 

ffoo = -/ + ^srn 2 0>O. (22) 

We have studied the appearance of ergoregions for a = 0.1, u>\ = uj2 and two different values 
of A, namely A = and A = 1.5, respectively. Our results are shown in Figl6j Very similar to 
the case for fixed uo\ = L02 an d varying a, the spherically symmetric, non-rotating boson star 
disappears from the solution for u\ < uJi iCr . We demonstrate this by plotting ^i(0) (4>i(r) has 
its maximal value at r = 0) as function of oj\. We find that 0J\cr decreases with increasing A, 
e.g. we find that u)\ iCr ~ 0.37 for A = and wi jCr ~ 0.35 for A = 1.5. This is easy to understand, 
since an increasing A means an increasing direct interaction between the two boson stars. We 
also plot the maximal value of 500 > 500,™ as function of uj\. Interestingly, for A = <?oo,m stays 
negative as long as the spherically symmetric boson star is present. Thus for a = 0.1 and A = 0, 
there are no instabilities due to ergoregions for those values of u\ = L02 where genuine interacting 
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boson stars exist. For smaller values of uj\ = to 2 (and the second branch of solutions), the curve 
for goo,m would follow the curve for a single rotating boson star given in [11]. The situation 
changes for A = 1.5. Here, goo,m becomes zero at a value of uj\ = which is larger that the 
value u>i :Cr at which the spherical boson star disappears from the system. We find u>i ~ 0.38 
for k = 0.1 and A = 1.5. Thus for toi^r < oj\ < to^ the two interacting boson stars possess 
an ergoregion signalling an instability. It appears that the stronger direct interaction between 
the boson stars tends to destabilize the system. A more detailed investigation and plots of the 
ergoregions will be presented in a future publication. 



4.1.3 h = 0, k x = 0, l 2 = 2, k 2 = 1 

This corresponds to the interaction of a non-rotating boson star and a rotating, axially sym- 
metric and parity odd - or in our notation "angularly excited" - solution. The (j) 2 component 
corresponding to the rotating boson star is odd under the z — > —z reflexion Our numerical 
results indicate that this case is qualitatively similar to the previous case. For vanishing and 
small a the energy density has its maximum in a ball centered around the origin as well as in 
two tori at z > and z < 0. This is clearly seen in Figd where we plot the energy density 
of the solutions for two different values of a, A = and uj\ = uj 2 = 0.8. When the parameter 
a increases, the spherically symmetric boson star disappears and the solution becomes purely 
axially symmetric. This is indicated in FigJ7]for a = 0.2 and a = 0.05. Again for a = 0.05, the 
energy density has a big value at the origin which results from the non-rotating boson star. For 
a = 0.2 this value has decreased significantly and the energy density looks like that of an axially 
symmetric solution. We also present the contour plots of the metric and matter field functions 
in FigE]for a = 0.1, A = and u\ = oj 2 = 0.8. 

Considering the ergoregions of these solutions we find that for a = 0.1 and both values 
of A = and A = 1.5, respectively, goo, m stays negative for all values of uj\ = ijj 2 for which 
genuine interacting boson stars exist. When interacting with a rotating, parity odd boson star 
the spherically symmetric boson star disappears from the solution for values of uj\ larger than in 
the case of interaction with a parity even boson star (see previous section). For uj\ = uj 2 smaller 
than wi jCr , i.e. the frequency at which 4>i(0) (the maximal value of <^>i(r)) vanishes (and for the 
second branch of solutions) the curve follows thus the curve of a rotating, parity odd boson star 
presented in |llj . We find that ui\^ C r ~ 0.43 for k = 0.1 and A = 0, while ui\^ C r « 0.41 for k = 0.1 
and A = 1.5. 
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5 Conclusion 



In this paper we have studied angularly excited as well as interacting Q-balls and their gravitating 
counterparts, boson stars. When a non-rotating boson star interacts with a rotating boson star 
for fixed angular frequency lo and varying gravitational coupling the non-rotating boson star 
tends to disappear from the system if the gravitational coupling reaches a critical value. For 
gravitational couplings above these critical values, only single, axially symmetric boson stars 
exist. The same holds true if the gravitational coupling is fixed and the angular frequency 
lowered. Then below a critical value of the angular frequency u the non-rotating boson stars 
has disappeared from the bound system. As a consequence the bound state of two boson stars 
with l± = ki = and 1% = k% = 1 exists only on a finite domain of the parameter space ui, UJ2- 
We observe that this behaviour is qualitatively independent of the choice of the direct interaction 
parameter A in the potential. 

When considering the existence of ergoregions for our solutions, which would signal an in- 
stability, we observe that the two cases of a non-rotating boson star interacting with a rotating, 
parity even boson star on the one hand and with a rotating, parity odd boson star on the other 
hand, are qualitatively different. While in the latter case no ergoregions appear for genuinely 
interacting boson stars, i.e. with the non-rotating boson star present in the system, ergoregions 
appear in the former case if the interaction parameter A is made large enough. 
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Figure 3: The contour plots of the metric functions f, I, g and m as well as of the scalar field 
functions 4>i and <p2 are shown for two interacting boson stars with l\ = 0, k\ = (spherically 
symmetric, non-rotating) and I2 = 1, &2 = 1 (axially symmetric, rotating, parity even). Here 
a = 0.2, A = and u\ = UJ2 = 0.8. Note that we use cylindrical coordinates z = rcosO and 
p = rsinO. 14 
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Figure 4: The dependence of the mass M, charge Q and </>i(0) on a is shown for two interacting 
boson stars with l\ = 0, k\ = (spherically symmetric, non-rotating) and Z2 = 1, &2 = 1 (axially 
symmetric, rotating, parity even) (solid). Here A = 0. For comparison we also give the values 
for the l\ = 0, k\ = single boson star for which = (dashed). 
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Figure 5: The energy density of two interacting boson stars with l± = 0, k± = (spherically 
symmetric, non-rotating) and I2 = 1, ki = 1 (axially symmetric, rotating, parity even) for 
a = 0.1 (upper curve) and a = 0.2 (lower curve), respectively. Here A = and u\ = UJ2 = 0.8. 
Note that we use cylindrical coordinates z = rcosO and p = rsin6*. 
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Figure 6: The value of the scalar field function at the origin associated to the spherically 
symmetric boson star ^i(O) as well as the maximal value of the time component of the metric 
<7oo,m are shown in dependence on uj\ = UJ2 for a = 0.1 and A = (solid) and A = 1.5 (dashed), 
respectively. Here l\ = 0, k\ = 0, h = 1, &2 = 1- 



17 




Figure 7: The energy density of two interacting boson stars with l± = 0, k± = (spherically 
symmetric, non-rotating) and li = 2, &2 = 1 (axially symmetric, rotating, parity odd). Here 
A = 0, uj\ = u)2 = 0.8 and a = 0.05 (top) and a = 0.2 (bottom), respectively. Note that we use 
cylindrical coordinates z = r cos 6 and p = r sin 9. 
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Figure 8: The contour plots of the metric functions /, I, g and m as well as of the scalar field 
functions (f>\ and <p2 are shown for two interacting boson stars with l\ = 0, k\ = (spherically 
symmetric, non-rotating) and I2 = 2 and A>2 = 1 (axially symmetric, rotating, parity odd). Here 
q = 0.1, A = and oj\ = L02 = 0.8. 
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Figure 9: The value of the scalar field function at the origin associated to the spherically 
symmetric boson star ^i(O) as well as the maximal value of the time component of the metric 
<7oo,m are shown in dependence on uj\ = UJ2 for a = 0.1 and A = (solid) and A = 1.5 (dashed), 
respectively. Here l\ = 0, k\ = 0, h = 2, k2 = 1. 
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